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Abstract 

We study the interfaces' time evolution in one-dimensional bistable extended dy- 
namical systems with discrete time. The dynamics is governed by the competition 
between a local piece-wise affine bistable mapping and any couplings given by the 
convolution with a function of bounded variation. We prove the existence of travel- 
ling wave interfaces, namely fronts, and the uniqueness of the corresponding selected 
velocity and shape. This selected velocity is shown to be the propagating velocity 
for any interface, to depend continuously on the couplings and to increase with 
the symmetry parameter of the local nonlinearity. We apply the results to several 
examples including discrete and continuous couplings, and the planar fronts' dy- 
namics in multi-dimensional Coupled Map Lattices. We eventually emphasize on 
the extension to other kinds of fronts and to a more general class of bistable ex- 
tended mappings for which the couplings are allowed to be nonlinear and the local 
map to be smooth. 
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1 Introduction 

The fronts in bistable discrete dynamical systems defined on a one-dimensional lat- 
tice have been investigated in various models such as in Lattice Dynamical Systems 
(LDS) |@], or in Coupled Map Lattices (CML) j| §. These models can be viewed 
as being in the class of (extended) dynamical systems for which the time evolution 
is given by 

u t+1 = Cu l + £ F o u\ (1) 

where u is a lattice configuration at time i e N (or Z), and C and £ are some linear, 
continuous and homogeneous operators, the so-called couplings. The mapping F is 
the direct product of a (local) bistable one-dimensional map /, i.e. the dynamics 
induced by / has two stable fixed points. In LDS C is often chosen to be the 
Laplacian to model a diffusive process and £ is a multiple of the identity ||, [ljj . 
On the opposite, in CML the diffusion is included in £ and £ is generally assigned 
to zero [H ||, H |j| . In all these models, the dynamics is governed by the competition 
between the fixed points along the lattice and then results in pattern formation. 
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Figure 1: The local map /. 



To describe the fronts in these space-time discrete systems, one has to consider 
bounded configurations depending on a spatio-temporal continuous variable, the 
so-called (front) shapes (5). It is then useful to extend the domain of L and of CI , 
and therefore the phase space of these extended systems, to bounded functions of a 
real variable. By doing so, one obtains a more general class of extended dynamical 
systems with discrete time for which, to each point of the real line is associated a 
bistable map, and these maps interact through C and This new class contains 
in particular LDS and CML. 

In this work, to obtain explicit expressions, the local map is chosen to be piece- 
wise affinc. However we shall see how the existence of fronts can be stated for 
some bistable models with a smooth, say C°°, local map. Up to a linear change of 
variable, it is always possible to assume the two local fixed points to be and 1. 
The resulting expression is then (see Figure 1) 

f(x)=ax+(l-a)H c {x), (2) 

where a £ [0, 1), c € (0, 1] and 



H c {x) 



if x < c 

1 if x 5s c 



is the Heaviside function at c.[] 

By combining relations (|l|) and (^|), one obtains the general expression for a 
piece-wise affine bistable extended mapping 

u t+1 = L lU l + L 2 H C o u\ (3) 

The goal of this article is to describe the dynamics of fronts and of interfaces in 
this system. Some general conditions are imposed (below) to L\ and to Li in order 
to ensure the model to be interpreted as a bistable mapping with homogeneous 
couplings. We shall also see how the results apply to the description of planar 
fronts in multi-dimensional CML. 

The plan of the article is as follows. We begin by defining the dynamical system 
and particularly the couplings under consideration. Still in Section 2, a represen- 
tation of couplings in terms of convolutions is given and several related properties 

1 In all the paper, stands for the Heaviside function at ui. 
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are reminded or established. The definition and the existence of fronts are exposed 
in Section ||. In this section, we also investigate the set of possible fronts veloc- 
ities and the set of the symmetry parameter c for the (non-) existence of fronts. 
Furthermore by introducing a distance in the set of functions associated with the 
couplings, the selected velocity is shown to vary continuously with the latter. In 
the following section, we prove that any interface, i.e. any orbit for which the state 
at each time is a configuration linking two different phases, always propagates with 
the fronts velocity. Section 5 deals with applying this formalism to the analysis of 
the fronts velocity in a LDS and in the continuous purely diffusive model. Moreover 
we describe how the planar fronts' dynamics in multi-dimensional CML enters in 
this framework. Finally the results stated in the plateaus, i.e. the intervals in c for 
which the velocity is constant, are extended to bistable extended mappings with 
nonlinear couplings and smooth local maps. We also prove the existence of other 
fronts as the unstable ones. 



2 Definitions and preliminary results 

2.1 The dynamical system 

First of all, the phase space is chosen to be the space, called Ai, of bounded Borel 
measurable functions on R, endowed with the uniform norm 

||tt[| = sup 

In this phase space, the action of the dynamics on the orbits {w'jt^o is given by the 
relation (||) where L\ and Li are linear continuous operators satisfying the following 
conditions 

(i) L\ and Li are s-homogeneous 

(ii) L\ and Li are positive 

(Hi) (Li + i 2 )l = 1 and L 2 + 0, 

where 1 stands for the function u(x) — 1. For wgR and u G Ai, let 

(a u u) (x) = u(x -to) Vx G R, 

be the shift (or translation) operator. An operator L with domain Ai is said to be 
homogeneous if 

La" = a"L Vcj G R, (4) 

and is said to be s-homogeneous if it is homogeneous and commutes with the point- 
wise limit of equi-bounded sequences in Ai, i.e. if {u n } n <z®, u n G Ai is such that 

3m G R : sup ||u n || < fn and Vx G R lim u n (x) — u(x), 

then 

Vx G R lim (Lu n ) (x) = (Lu) (x). (5) 

n — > + oo 

The s-homogeneity reflects the homogeneous and the punctual characteristics of 
the dynamics An operator L is positive if Vw £ Ai such that u 0, we have 
Lu ^ 0.0 The condition (ii), together with (Hi), indicates that the dynamics (||) can 
be interpreted as the competition between a local reaction and a global diffusive 
process. The condition (Hi) guarantees the invariance of the local fixed points by 
the couplings. 

2 u ^ means u(x) ^ Vx e R. 
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After having denned the dynamics using s-homogeneous operators, we now 
obtain an integral representation of the latter using the Lebesgue-Stieltjes inte- 
gral. The following statement reveals the existence of a bijection between the 
s-homogeneous operators on M and the right continuous functions of bounded 
variation vanishing at — oo. This set of functions will be denoted by BVq. 

Proposition 2.1 A linear and continuous operator L of M. into itself is s-homogeneous 
iff it exists a unique function h G BVo, such that 

VuEM, {Lu){x) = J u(x- y)dh{y) Vx G M. (6) 

Moreover L is positive iff h is increasing. 

Proof: If the integral (^) defines L, then the Lebesgue dominated convergence 
theorem shows that L is s-homogeneous. Conversely let h(x) = LHq(x). Then by 
relations (H) and (nh, h is right continuous and lim h{x) = 0. By using (El) and 

x — > — oo 

the linear continuity of L, we conclude that h is of bounded variation in K (similarly 
as in the proof of the Riesz theorem in jlO]]). Then (^J) is true for all the functions 
of the form 

N 

u=/ J s n (H Xn _ 1 - H Xn ) . 

Any continuous function is the point-wise limit of a sequence of functions of this 
form. Relation (||) then holds in the (Baire) class B° of all bounded continuous 
functions. By the same argument, it is also valid in the class B 1 of functions which 
are point-wise limits of functions in B°. We can proceed in this way transfinitely to 
conclude that the proposition is true for any bounded function in a Baire class B a 
(where a is a countable ordinal). But 1J B a is exactly the set of all Borel measurable 

a 

functions Q.| □ 
Let us recall that, for a function h £ BVo, the convolution is defined by 

Vw G M h* u(x) = J u(x - y)dh(y) Vx G R. 

We have adopted an unusual ordering in this definition to indicate that the con- 
volution by h G BVo corresponds to the action of an operator in M. Moreover 
notice that the convolution h * it is defined even if u ^ BVo and is thus in general 
non-commutative. 

We then obtain an equivalent expression for the dynamical system (||) 

u t+1 = h 1 *u t + h 2 *H c ou t , (7) 



where hi (resp. h 2 ) corresponds to L\ (resp. L 2 ) according to Proposition 2.1 and 
is thus an increasing function in BVo. We denote by T the subset of BVq and of 
M. composed with right continuous increasing functions vanishing at — oo. This 
alternative expression for the dynamical system, combined with the properties of 
the convolution and of the functions in I 1 both described in the following section, 
will prove to be convenient for the analysis of fronts. To conclude this section, 
we mention that the interfaces' dynamics in a similar but space-time continuous 
bistable extended dynamical system, i.e. a first order in time differential equation 
for which the interaction is given by the convolution with a C 1 function, has been 
investigated in ||. 

3 Requiring the functions to be bounded is not restrictive because any unbounded Borel measur- 
able function is a point-wise limit of bounded Borel measurable functions. Hence we can consider 
the classes B a containing only bounded functions. 
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2.2 Basic properties 

Using Fubini theorem, it is shown that the convolution is commutative and asso- 
ciative in BVq. As a consequence, two any linear continuous and s-homogeneous 
operators commute. 

The following properties are a direct consequence of the definition 

(i) If h\,h% G X, then hi * h<z £ X and \\hi * h<z \\ = \\hi || ||/i2||- 

(ii) Let fcel and ui,M2 £ Ai. If u\ ^ then h * ui ^. h*U2- 
(Hi) \/u E M cr^u = *u, 

(iv) Vw, w' € K ff w * = • 

Moreover let the left closure of a real subset A be defined by0 

A £ = f| A+[0,S). 

<5>0 

Using the diagonal process, one shows two basic properties of the left closure. 
1 

I l 

(v) A +B =A + B . 

(vi) ~\JJ?*= \JA*. 

We now introduce two sets which (partly) characterize the right continuous 
increasing functions. For the interpretation in terms of coupling, they determine 
the sets of points coupled to and thus, using the homogeneity, the generalized 
lattice on which the coupling acts. For ftcl the set of increase points is defined 
as follows. 

E(h) = {x e K : V5 > h{x) > h(x - 5)} . 
The set of discontinuity points, a subset of E(h), is given by 

D(h) = {x e M. : h(x) > h(x - 0)} , 

where h(x — 0) — lim h(x — S). As a first result, one has a sufficient condition 

5-+0,5>0 

for E(h) to reduce to D(h). 

Lemma 2.2 and E(h) is countable then E(h) = D(h). 

Proof: If ft. G I, then it is equal to the sum of a continuous function h c and a step 
function hd which are such that h c , hd S X and D(h) = E(hd) = D(hd) |lC{ ]. Hence 
E(h) \ D(h) C E(h c ). Now by assumption let {a„}„ e N be an enumeration of E(h c ). 
We obtain 

/ dh c = ^2 h c {a n ) - h c (a n - 0) = 0, 

because Hh{E{h) c ) = 0, where \Xh is the Lebesgue-Stieltjes measure associated with 
h <E X, and where A c denotes the complement of A in R. As a consequence h c = 
and h — hd- □ 
Moreover, the evolution of these two sets under the convolution and the sum- 
mation is given in the following statement. 

Proposition 2.3 (1) Let h, h' G X. Then h * h' G X and 

(la) E{h*h') = E(h) + E(h') e , 
(lb) D(h*h') = D(h) + D(h'). 

4 i.e. if ui(x) ^ 1*2(2;) Vx 6 K 

5 A + B = {x + y : x e A and y £ B} 
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(2) Let {/ii}i g N be a uniformly summable sequence of functions in I. Then h 
/^hj G T- and 



(2a) E(h) = \jE(hi) 
(2b) D(h) = \jD(h t ). 



Proof: We only show relation (la) here. Relation (lb) can be proved similarly. 
Statement (2b) is a direct consequence of uniform convergence and (2a) follows 
from (|8|) below. 

For an increasing function, the definition of E{h) implies that for a given 8 > 0, 

h{x) -h{x-S)>0 iff (x -S,x]n E{h) ^ 0. (8) 
Consequently given h,h' Gl we have for any x € R and 6 > 

h * h'(x) — h * h'{x — 5) = I h(x — y) — h(x — y — 5)dh'(y), 

J {-E(h)+(x-8,x])nE(h') 

from which we deduce the relation 

\/6>Q h*h'(x) -h*h'(x-S) > iff V5 > x e E{h) + E{h!) + [0, 6). 
This equivalence is nothing but (la). □ 

3 The fronts 

With the previous definitions and properties provided, we are now able to analyse 
the dynamics of the simplest structures linking two different phases^], namely the 
fronts. In this section the latter are firstly defined. Then we state their existence, 
in particular by considering various conditions on the couplings hi and hi and 
on the parameter c. We also prove the continuous dependence of their velocity 
on the couplings using an adequate topology. We finally consider the existence of 
anti-fronts by relating it to the fronts' existence. 

3.1 The existence of fronts 

Definition 3.1 A front of velocity v is an orbit of (j?|) defined by 

u* = H vt+X0 *<f> yt G z, 

where Xq £ M and the shape <fi <E M. is such that 

H c o0 = H Q . (9) 

This definition is somehow restrictive because, for the sake of simplicity in this 
section, it imposes 0(0) ^ c. However, we shall see in Section ^ how to relax this 
condition by imposing the shape to satisfy, instead of (|^), the condition cj)(x) < c if 
x < and <j){x) ^ c if x > 0, and therefore by letting 0(0) arbitrary. 

Because of the translational invariance in (Q), if it exists a front, then there is 
also an infinite number of fronts with the same velocity and the same shape, one 
for each value of xo . 

6 The phases of the system are the fixed points and 1 which compete together when coupled. 
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As we shall see below in the proof of Theorem 3.2, this definition implies that 
the shape <f> is a right continuous increasing function which depends on the velocity. 
Moreover it behaves asymptotically as 

lim 4>{z) = and lim cf>(z) = 1. 

z — > — oo z — >+oo 

The fronts thus defined are actually the travelling waves resulting from the motion 
of an interface between the two different phases of the dynamical system under 
consideration. 

The natural questions about the fronts deal with their existence, their uniqueness 
and the selected velocity given two any couplings hi and hi and a value of the 
symmetry parameter c. It is also of interest to determine the set P of possible 
velocities when c varies in (0, 1], and the couplings are kept fixed. 

Theorem 3.2 Given two couplings hi and hi, there is a Baire first category set G 
such that for any c G (0, 1) \ G, the dynamical system has a unique front shape. 
The corresponding fronts have a velocity v(c) which increases with c and which can 
be extended to a left continuous function on (0,1). The set P of possible velocities 
is a dense subset of the range of the function v and has the following bounds 

inf P = v min and sup P = w max , 

where v m - ln — min {inf E(hi), inf E(hi)} and w max = max {sup E(h±), sup EQ12)}. 

Assuming some additional (general) conditions on hi and hi, the set G is shown 



to be countable and nowhere dense (see section 3.4). If v m i n = inf E (hi) and 
» ma x = sup E (hi), then P is dense in (v m i n! Vmax) and therefore v is a continuous 



function of c (see Section 3.3). If the function hi is continuous, then G is empty, P 
contains (v m in, t) max ), v is continuous and the shape cj> also varies continuously with 
its argument. 

For c G G no front exists. We have instead a "ghost front", i.e. a sequence of 
configurations in T which is not an orbit of (|7]) but which attracts a set of initial 
conditions || . Ghost fronts are due to the discontinuous character of the local map 
possibly combined with the couplings. The ghost front's velocity v(c) may belong 
to P but it may not since, in general, P and the range of v(.) are different. If v(c) 
belongs to P this ghost orbit can be avoided by changing the value of the local map 
at c for fixed couplings. 

3.2 Proof of the fronts' existence 



The proof of Theorem 3.2 consists in investigating the function <j>, which can be 
defined independently of the front existence, to deduce the conditions for it to be a 
front shape. 

Assuming the existence of a front of velocity v, we obtain from (Q) 

H v *(j> — hi*(p + hi* Ho. 

Then by using the convolution properties, we solve this equation to obtain the 
following expression, where the dependence on v has been added 

00 

(/>(z,v) =^h 2 *h*i k (z+(k+l)v) VzgM, (10) 

fe=0 

where for any h G BVq we denote 

h*° = H and h* [k+1) = h* k *h for k > 0. 



7 



Consequently a front of velocity v exists iff the function defined by ([H]) satisfies 
condition (^). From relation (10), the convolution properties, the condition (Hi) on 
the couplings and Proposition 2.3, it turns out that 

V»eK 4>{.,v)el, (li) 

and 

1 

oo 

E(<f>(.,v)) = (J E{h 2 ) + + k E{h 1 ) + {-(k + l)v} , 

where for any real subset A we use the notation 

+ A = {0} and +(fc+i) A = + k A + A for k ^ 0. 
Similarly we have 

Vzel (12) 

and 

- k E{h x ) 

fc + 

k=Q 

where for any j/£K\ {0} 



— = {- : x e A}. 

y y 

Property ([ll]) gives the following necessary condition for a front of velocity v to 
exist for some c S (0, 1] 

G E(cj>(.,v)). (13) 

Actually if £ E(4>(.,v)), then 

35 > : Ve G (0,(5) 0(-e, «) = 0(0, «), 

and the condition (j^) cannot be satisfied for any c. Conversely by choosing c = 
0(0, v), we conclude that for all v for which relation ( fl3| ) holds, there is (at least) 
a value of c for which the fronts of velocity v exist. The set of possible velocities is 
then given by 

P = {veR : e E((j)(., v))} . 
This set has the following properties 

U £ % + + T (M CPC W (»,.», (14) 

fe=0 

from which the bounds on P are deduced.^] 

Now one has to determine which fronts of possible velocity are selected given 
a value of the parameter c. From relation (10) it turns out that if v\ < i>2 then 
0(0, v±) ^ 4>{—{v2 — vi),V2), and thus if c ^ 0(0, Ui) then the existence condition 
(||) shows that V2 cannot be the fronts velocity for this value of c. According to this 
property, the velocity of the (assumed) existing fronts is uniquely given by[] 

u(c)=min{ueR : 0(0, v) ^ c} . (15) 



7 Sinco E(h) = E(h) , we have sup E(h) £ E{h) if sup E{h) ^ oo. Consequently u max 6 P iff 

) — su 

fc+1 



s\ipE(h\) ?C sup_E(/i2) < +oo and in all (finite) cases v = supE(hi) + sup fi ( fe 2> sup g(fci) ^ p 



for any k ^ 0. 

8 i}(c) is defined at c = 1 iff v ma x < oo. Indeed, if t) max = oo, the set (v £ M : 0(0, u) l| is 
empty. 
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The corresponding fronts exist iff 



V6>0 (j>(-S,v(c)) < c, 
and the non-existence set is given by 

G = {c : 35 >0 c = (f)(-8,v(c))}. 

From its definition and from relation (^) the function v is clearly increasing. To 
show that it is left continuous, we notice that given S > one has (f>(0, v(c) — S) < c. 
Let then c' be such that 0(0, v(c) — 5) < c' < c. Again from the definition, we obtain 
v(c) — S < v(c') and consequently v is left continuous since it is increasing. Moreover 



its range is E((j)(0, .)) and thus relation (14) shows that P is a dense subset of this 
range. 

To achieve the proof we have to show that G is of Baire first category. To that 
goal let 

G (n) = (c : c = 6(---0,v(c)) 
(_ n 

We have 

oo 

G= \J G {n \ 

71=1 

and it then remains to show that each is nowhere dense. By contradiction, we 
assume that G*-™* 1 is dense in the (non-degenerated) interval (ci, C2). Then (ci, c%] C 
G^ n ' since both the functions v and 4>(—— — 0, .) are left continuous and increasing. 
It follows that (ci, C2] C G. Now for any c £ (ci, 02) we have 

0(0 - 0, v( Cl )) < ci < c = 0(0 - 0, «(c)) < c 2 = 0(0 - 0, v(c 2 )), 

that is to say v(c) € (w(ci), ^(02)) and «(c) ^ P. This is impossible since v is 
increasing and P is a dense subset of the range of v. □ 

3.3 Characterization of the velocity set P 

We now proceed to the analysis of P and to its consequences. In the following 
lemma we firstly describe P inside and outside the convex hull of E(hi). Let us 
define 

ij = inf E(hj) and Sj = sup E(hj) for j = 1, 2. 



Lemma 3.3 (i) If #E(h\) > 1 then for any 9 £ (0,1), 6 irrational, we have 

Va,beE(hi), a + 9(b-a)eP. 
(ii) If S\ < S2 < +00 then 

Pn( Sl ,s 2 ] =-E(0(o,.))n(si,s 2 ]. 

(Hi) If #E(hi) > 1 and s 2 = : +00 t/ien for any 9 > 1, 9 irrational, we have 
Va, b £ £(/ii), a < 6 : a + 0(6 - a) e P. 
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The proof is given in Appendix |A|. Naturally one has statements similar to (ii) 
when — oo < 12 < i\ and to (Hi) when 12 — —00. Namely 

(ii') If —00 < %2 < i\ then 

Pn[i 2 ,* 1 )=B(0(o,.))n[*2,ii). 

(Hi') If #E(hi) > 1 and 12 = —00 then for any 9 > 1, 9 irrational, we have 
Va, 6 G a<b : b + 9{a - b) G P. 

Moreover since v is a left continuous increasing function, a discontinuity of v 
means the absence of an interval in its range. As P is contained in this range, 
statement (i) of the previous lemma implies the following. 

Corollary 3.4 In the interval u _1 ((ii,si)), the selected velocity v(.) is a continu- 
ous function of c. 



Furthermore still in statement (i) of Lemma we have seen that all the 
irrational velocities in (i\, s\) are fronts velocities. We now give a sufficient condition 
for all the rational velocities in this interval to be fronts velocities. To that goal the 
lattice generated by a real subset A is introduced 



R(A) = A+ \J+ k (A-A), 

k=0 

where A denotes the (usual) closure of A. Some characterizations of the lattices 
generated by real subsets are listed in the following remark. 

Remark 3.5 (i) If R(A) / R and f^A > 1, then it exists a G M + such that 
R(A) = a + ouL for any a G A. 

(ii) If A contains an infinite bounded set, then R(A) = R. 
(Hi) If there are x , xi,x 2 G A such that x x \Z x x ° a G M \ Q, then R{A) = R. 
(iv) IfE(h) ^ D(h), then R(E(h)) = R (see LemmaK_ 



Proposition 3.6 If E(h 2 ) D R(E(hi)) ^ 0, then («i,si) C P. 



In particular, if one of the last conditions in Remark 3.5 holds for E{h\) then 
(ii, si) C P independently of /12. 

Proof: Let d G E(h 2 ) (~l R(E(hi)) and 5 > 0. There are I ^ 1, {a„} ^n^z, a n G 
E(h±), a n < a n+ i, and {mi}o^i<ij mi G Z, such that 

l-i ^ 
d a + ^TOj(a i+ i - a,i) < d+ -, (16) 

2 = 

We can suppose that oq is arbitrarily close to i\ and a/ = Si by choosing, if necessary, 
too = 1 ai 

such that 



too = 1 and mi_i = 0. Now let J = 2 max |m»| and if I > 1 let M.N G Z + be 



i-2 j. 





< N{ ai - o,_i) - M ^(Z - i) J(a i+1 -<»*)<-. (17) 

i=0 

If / = 1 we set N = M = 0. Let us choose any < p < q co-prime integers and 
define 

ni-i = —mi-i — N + jp and = -m; + M(l — i)J + jp for ^ i ^ I — 2, 
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where the integer j is chosen large enough so that 

jq — 1 ^ n ^ n x ^ • • • > nj_i > 0. 
By adding © and (0) we obtain 

rf«:«o- - 0+ l)-)(a i+ i - a,) < d + 6, (18) 

Q 

where k = jq — 1. It is easy to see that the latter is equivalent to 

-S < d+ (k - n Q )a + (n - ni)aH h (rt;_ 2 - n/_i)a/_i + n/_ia/ - (k+ l)v < 0, 

where u = ao + - (a; — ao). Since <5 is arbitrary, we deduce that ugP and the proof 



is achieved by using Lemma 3.3 (i). □ 



3.4 Characterization of the non-existence set G 



Phrased differently, Theorem 3.2 claims to each velocity v G P corresponds a set 
of the parameter c for which the only existing fronts are those of velocity v(c) and 
their shape is unique, the sets corresponding to different velocities being disjoint. 
We now describe both these sets and the complement of their union, i.e. the set G. 

Firstly if v(c) e P \ D(<J)(0, .)) then 

4>{x, v) < cj)(0, v) = c Vx < 0, 

and the fronts exist. If v(c) £ D(cj)(Q, .)) since we have 

0(0 -0,v) = 0(O,v-O), 

then either 0(0 — 0,v) < c and the fronts exist, or 0(0 — 0,v) — c. In the latter case, 
condition (M) is satisfied or not according to v(c) e Pi or v(c) ^ Pi where 



Pl = {«eM : 06fi(^(,(j))\{0}'} 



In this way, we conclude that for h x , /12 and v G P fixed, there exists a set of the 
parameter c given by 

c=0(O,u) if v G P\L»(0(O,.)) 

0(0, v - 0) < c < 0(0, «) if u G r>(0(O, .)) (19) 
c = 0(O,u-O) if v G D(ct>{0, .)) n Pi 

for which (only) the fronts of velocity w(c) = w exist. 

In particular, if for a given velocity v G P, it exists an interval such that the 
selected velocity is v(c) = v for c in this interval, w is said to correspond to a 
plateau. According to the previous conclusion, a velocity v corresponds to a plateau 
iff v G £>(0(O, .)). Notice that Proposition applied to relation ( p~p| ) leads to 



D(0(O,.))= U 



D(h 2 ) + + k D{hi 



k + l 

k=0 

Let us now give some general conditions for the couplings to satisfy in order 
for G (the set of c values for which no front exist) to be countable and/or nowhere 
dense. According to the previous reasoning, this set is given by G = Gi U G2 where 

G 1 = {0(O,t>) : «GP(0(O,.))\P} and G 2 = {0(0 - 0, v) : v G Z)(0(O, .)) \ Pi} ■ 

G2 is countable since the discontinuity set of an increasing function is ppj . It 
is nowhere dense by the inequalities (|l^). Hence it remains to investigate these 
properties for G\. 
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Proposition 3.7 G is countable if one of the following conditions holds 

(i) #E(h{) ± 1 

(ii) #E(hi) = 1 and E(h 2 ) is bounded. 

It is possible to find some examples in which G is not countable. 

Proof: Firstly if hi = then (j>(x,v) — h 2 (x + v). In this case we obviously have 
E(4>(0, .)) = P and thus Gi = 0. Assume now that #E(hi) = 1 and E(h 2 ) is 
bounded. Then statement (ii) of Lemma |3.3| implies that #Gi ^ 1. Finally if 
#E(hi) > 1 then (ij.,Si) \ P is countable as follows from statement (i) of Lemma 



3.2. So is (E((j)(0, .)) \P)\ [ii,si] as it can be deduced from the assertions (ii) and 
(Hi) of the same lemma. □ 

To state conditions under which G is nowhere dense in (0,1), it is useful to 

introduce a notation and some definitions. Given two real subsets A and B, the 

00 

notation A C B has the following meaning 



AcB iff 



sup A = +00 sup B = +00 
inf A = —00 =>■ inf i? = —00 



In particular, if A is bounded no restriction is imposed on B and thus A C B for 
any Bel. The second concept that will be of use is the following. 

Definition 3.8 A subset A o/R is commensurable if 

X\ — Xq 

Vxq,xi,x 2 €A x 2 f= x =>• € <y). 

x 2 - x 

A is said to be incommensurable otherwise. 



Proposition 3.9 G is nowhere dense if one of the following conditions is satisfied 

(i) E(hi) is incommensurable 

(ii) hi = 

OO 

(Hi) #E{hx) = 1 and E(h 2 ) C D{h 2 ) 

OO 

(w) E(h 2 ) n R(E(hi)) ^ and E{h 2 ) C D(h 2 ) U E^) 
(v) D(h 2 ) j= and E(h 2 ) C D(h 2 ) U S(fti). 

Although many cases do not fit in these propositions, we can conclude that G is 
countable and nowhere dense in most of the physical examples of piece-wise affinc 



bistable extended mappings. Actually, Propositions 3.7 and |3.9| cover the case of 
cellular automata (i.e. hi — 0), the case of lattice dynamics (i.e. D(h 2 ) = E(h 2 )) 
(with finite range for G to be countable) from which belong LDS and CML, and 
more generally the case E{h 2 ) C E(hi) that is expected to occur in the model ([!]) 
since we have 

£u + C'F ou= (£ + aC) u + (1 - a)C'H c o u, 
and thus C + aC is expected to act at all the points where £ acts. 



Proof of Proposition |3.fi| : (i) If E{h\) is incommensurable then Proposition p.6| 
implies (ii, si) C P. Now if s 2 = +00 let xq,Xi,x 2 E E(hi) be such that ^Z^° E 
K\ Q. By using statement (Hi) of Lemma |3.3| with a — xq, b — Xi and with a = Xq, 
b = x 2 respectively, we conclude that (si, +00) C P. If s 2 < +00 then Lemma |3.3| 
(ii) results in (si,s 2 ) R E(<j)(0, .)) \ P = 0. One can proceed similarly for the cases 
i 2 = —00 and i 2 > -co to conclude that E(<p{0, .)) \ P has at most two points (ii 
and si) and consequently #Gi ^ 2. 
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If h% = then Gi = (see the proof of Proposition 3.7) and (ii) follows. 

Since we have already proved (i) and (ii), we can now assume that E(h\) is com- 
mensurable and non-empty. Hence it is countable and using Lemma |2.2| , we have 
E(hi) = D(h±). The proof of the remaining assertions is based on the following 
lemma proved in Appendix |^. 

Lemma 3.10 (i) If D(h 2 ) 7^ and E(hi) = D(hi), then G is nowhere dense in 
V _1 ((*l> s l))- 

(ii) If sup D(h 2 ) = +00 and a G D(h\), then G is nowhere dense in v~ 1 ((a, +00)). 
(Hi) //inf D(h 2 ) = —00 and a G D{h\), then G is nowhere dense in v _1 ((— 00, a)). 



To achieve the proof of Proposition 3.E, we are going to show that under the con- 
ditions in (in), (iv) and (v), G\ is nowhere dense in each of the following intervals 
V {(i 2 , *i)), v~ 1 ({ii,s{)) and v~ 1 ((si, s 2 )), concluding that G\ and consequently 
G is nowhere dense in (0, 1). 

00 

Firstly we prove that the condition E(h 2 ) C D[h2)LiE(hi) implies that G is nowhere 
dense in v _1 ((si, s 2 )) by considering the following cases. If supD(h 2 ) = +00, then 
we use Lemma |3.10] (ii). If supE(hi) = +00, then there is nothing to prove because 



v 1 {(si,s 2 )) = v 1 (0) = 0. If supE(h 2 ) < +00 Lemma 3.3 statement (ii) results 
in v- 1 ((s 1 ,s 2 ))nGi =0. 

Similarly G is nowhere dense in v^ 1 {{i 2 , 

Now if #E(hi) = 1 then w _1 ((zi,si)) is empty, proving (Hi). If E(h 2 ) r\R(E(hi)) ^ 



then by Proposition 3.6, v 1 ((ii,si)) n G± = and (iv) is proved. Finally if 



D(h 2 ) 7^ then by Lemma 3.10 (i), G is nowhere dense in v ((ii,si)) and ('uj 



holds. □ 



3.5 Continuity of the fronts velocity 

We now consider the dependence of the fronts velocity on the couplings. In par- 
ticular, we are going to show that small variations of the couplings induce small 
variations of the selected velocity. As couplings variations, we also would like to 
allow for changes in the sets of increase points of the couplings hi and h 2 . For in- 
stance, this is the case when the planar fronts' direction varies in multi-dimensional 



CML (see Section 5.3). We start with a definition of a distance in X. For h,h' G X 
let 

d(h, h') = inf {e > : h(x - e) - e < h'{x) < h(x + e) + e Vz G E} . 

One easily checks that d(., .) is a distance. For this distance, the ball of radius e 
centered at h is the set of functions for which the graphic lies in the band of width 
e in the direction of the line y — —x around the graphic of h. Similarly in the set 
of left continuous increasing functions defined on (0,1), we define d(v,v') as the 
infimum of the positive numbers e such that 

v(x — e) — e ^ v'(x) if x G (e, 1) 
v'(x) < v(x + e) + e if x G (0, 1 - e) ' 

From this definition one concludes that d(.,.) is a distance and d(.,.) ^ 1. The 
continuity of the selected velocity in the subsequent topology is an immediate con- 
sequence of the following statement. Let v (resp. v') be given by (|l5|) for the 
couplings h\ and h 2 (resp. h[ and h 2 ). 
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Proposition 3.11 For the couplings hi, h[, hi, h' 2 , let 

m = max{d(hi, h[), d(h2, h' 2 )}- 

Then 

d(v,v ) < 



^4s a consequence, for any fixed c <E Si)), v is a continuous function of the 

couplings in the Euclidean topology. 

Proof: Let us firstly show the continuous dependence of the shape in the topology 
induced by d(., .). We notice that for any e > d(h, h'), the following holds 

H e * h 1 - e ^ h < i?_ e * ft' + e. 

Hence by induction and by using the convolution properties, we obtain for any k ^ 1 

fe-i fe-i 
tf fce * - ll/illil/i'f" 1 -' < fe* fc ^ H_ ke * + H/ilhl/i'll*- 1 -*. 

i=0 i=0 

Replacing these inequalities in the expression of 4> and using again the convolution 
properties result in the following 

Ve > m 4>'{x, v-e)- - _ 2 ^, | <j>{x, v) < 0'(a;, v + e) + JZ^J ' 

where (f)' is computed with h\ and ftj. From the right inequality we obtain 

c -^W^ (0 ' 5(c)) -T^W^' (0 ' ,1(c, + e) ' 

and hence by the definition of v, it follows 

2e 



Fn 0(c) + e. 



Similarly we have 



c - yqj^ < ^(0, Ac)) - < 0(0, C(c) + e), 

and by the definition of v and by replacing c by c + 1 Jp l > | | , we obtain 

Hence d(v, v') ^ || anc ^ tne statement is obtained by choosing e arbitrarily 

close to to. 

Finally, we assume that v(c) e (ii, si). We have shown that 

2e \ 2e » . . / 2e \ 2e 



i> c 



i - K|| J i - 11*111 * c ' (c) * H c+ 1 - Kll J + i - Kll' 

for e > to. Thus from the continuity of v(.) in v~ 1 ((ii, si)), we obtain for c fixed 



lim v'(c) — v(c). 

m—>0 
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□ 



3.6 The anti-fronts 

To conclude the study of the existence of travelling wave interfaces in our dynamical 
system, we mention that the results about the fronts can be employed to deduce 
the existence of anti-fronts and the uniqueness of their shape. An anti-front is a 
travelling wave for which the shape cf> antl obeys the condition 

H c o (j> wti (x) = H (-x) \/x e R. 

By defining 4> a (x) = </> antI (— x) and by using the shape equation, one can easily 
prove the following statement. 

Proposition 3.12 The dynamical system (Qj has anti-fronts of velocity v and shape 
^anti £/j ere RX { s i fronts of velocity —v and shape (ft a for the dynamics given by 

u t+1 = ft** it* + h a * H c ou\ 

where hf(x) = \\hi\\ — hi(—x — 0) for i = 1, 2. 

Moreover it is possible to relate antl (c) the velocity of anti- fronts to the fronts 
velocity for the same couplings, i.e. those represented by hi and /12, in the following 
way 

anti (c) = 0(1 -c + 0). 

Hence in general the selected velocities of fronts and of anti- fronts, as well as the 
corresponding non-existence sets G and G antl , are different. 



4 Velocity of interfaces 

In this section, an also important result, when one has in mind physical applications 
of bistable extended mappings, is presented. In fact we are going to prove that the 
fronts velocity v is actually the propagation velocity for any orbit of (Q) composed 
at each time of a configuration linking two different phases, namely a c-interface. 

Definition 4.1 A c-interface is a function u G M. such that 

3 J- J e K ■ { U ^ <c if x <J[ 

" s [ u{x) ^ c if x > J s 

The set of c-interfaces is denoted by J c and Ji(tt) and J s (u) are the functions from 
J c to M defined by 

Jl(u) = inf{x : u(x) c} and J s (u) = sup{a; : u(x) < c}. 

The following proposition guarantees that the image under the dynamics (Q) of a 
c-interface is a c-interface. 

Proposition 4.2 If c S (0, 1) and u £ J c then h\ * u + hi * H c o u G J c . 

This statement is still true for c = 1 if v max < +00. 
Proof: Let u 1 = h\ * u + hi * H c o u and u be defined by 



u{x) 



u(x) if x ^ J s (it) 
c if x > J s (u) 
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We have u ^ u. By using the convolution properties, we then obtain 

liminf v}{x) ^ c\\hi\\ + \\h 2 \\ > c, 

x — >+oc 

since \\h\\\ = 1 — \\h2W and c < 1. Similarly, one proves 

limsupu 1 (x) ^ c||/ii|| < c. 

x^ — 00 

□ 

Hence the terminology c-interface orbit {w'} te N is justified. As it has been pointed 
out in the case of CML [||, these orbits can evolve towards different attractors. So a 
general common (asymptotic) stability result (in the sense of Lyapunov) in J c can- 
not hold. However the c-interface orbits all have a common dynamical characteristic 
which is the same propagation velocity. 

Theorem 4.3 If c is a point of continuity of the function v and {m'}( 6 n is a c- 
interface orbit, then its velocity exists and is given by v(c), that is to say 

hm = lim = v(c). 

t t^+ca t 

The result cannot hold if c is a discontinuity point of v because in this situation, 
the dynamical system has (also) ghost fronts of velocity v(c + 0) which attract some 
initial conditions. The proof of Theorem |4.3| is based on the following proposition. 



Proposition 4.4 For any fixed c G (0, 1) we have Vu° G J° Ve > 3(To,o"i G 
such that 

cr + t{v[c) - e) < Ji(u') s$ J s (u*) ^0-1+ tv(c + e). 

Clearly if u max < +00, we have a similar claim for c = 1. 

Proof: Given u° G J c \ by iterating the dynamics (0) we obtain for t ^ 1 



= *u° + J2 h ?* h *i k * H c o u t ~ 1 ~ k . 



fc=0 



Since {u*}t S N is composed of c-interfaces, the convolution properties and the defi- 
nitions of Ji and of J s lead to the following inequalities 

t-i t-i 

h\ i *u°+j2 h 2* h i k * #S.( U *-i-*) < «* < h *i * u ° + J2 hi * h * k * H uu^- k )^ 

k=0 k=0 

where H°(x) = H u (x - 0) Vw, ieE. Now given e > let 

A = c- 4>(0,v(c) -e) and A x = </>(0, v(c + e)) - c 

By the definition of w(c), we have Ao > and Ai ^ e > 0. Let then foj £i G N be 
such that 

||/n||*°||tt ||<Ao and ||/ ll || tl (|| W °ll+l)<A 1 , 
and (To, o~i G R be defined by 

(To = mm (Ji( uk ) ~ k(v(c) — e)) and o"i = max ( J s (u k ) — kv(c + e)) . 

We are going to prove by induction that 

(T +fc(u(c) - e) < Ji(u k ) < J s (w fe ) < cti +fcw(c + e) Vfc G N. (20) 
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By construction the left (resp. right) inequality is true for ^ k ^ to — 1 (resp. 
^ k ^ t\— 1). Assume the left (resp. right) inequality in (|2^) holds for ^ k ^ t— 1 
with t ^ to (resp. t^t%). It comes out that 



t-i 



u*(x) < |^i||*° + ^ ^2 * * i^o(^ - o-o - - 1 - fc)(«(c) - e)) V.t G R, 

fc=0 

and 

t-1 

u*(a;) ^ -\\h 1 \\ H \\u°\\ + ^h 2 *h* 1 k *H${x-<T 1 -(t-l-k)v(c+e)) Viel. 



fc=0 



Hence using the definition of the function <p(z,v) we have 

u\x) < ||/ii||* ||w°|| + <?!)(a;-cr -i(w(c)-e),u(c)-e) Vi e R, 
and for any S > 

u*(a;) ^ -H/iiH^dlw !! + 1) + 0(s - tri - «5 - t«(c + e), v(c + e)) Vx € R, 



fe=i. 



= || /nil* 1 . Conse- 



the latter being obtained by using the relation 
quently if x ^ <7o + t(v(c) — e) we have 

u*(x) < ||/ii|| to ||u ||+^0,S(c)-e) <c, 
in other terms do + t(,v(c) — e) ^ Jj(u*). If x ^ Ci + <5 + £v(c + e) then 

u\x) > -\\hi\\ H (\\u°\\ + 1) + 0(0, v(c + e)) > c, 
proving the right inequality J s (u*) ^ u\ + tv(c + e). u 



Proof of theorem 4.c: From the last proposition, it is deduced the following inequal 
ities 



Ji(u') Jb(V) 
Ve > w(c) — e ^ liminf — ^ limsup 

t-++oo £ 



and thus 



u(c + e), 



w(c) < liminf < limsup < w(c + 0). 



□ 



Since Jj(it*) ^ J s ( ut ) we obtain the existence of the limits if w(c) = v(c + 0). 

5 Examples 



In order to illustrate the results obtained in the general framework, we now present 
three examples. One focuses on discrete couplings, the second deals with continuous 
couplings, and we thirdly show how planar fronts in multi-dimensional CML enter 
in the general frameworkj^] 

9 In this section, we will refer both to the set P of poss ible velocities and to G the set of values 
of c for which no front exists. (See proof of Theorem 3.2) 
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5.1 One-dimensional Lattice Dynamical System 

In this first example, we choose the couplings as acting on a one-dimensional lattice 
and as being of finite range. The interest here resides in emphazing on two char- 
acteristics. One is the absence of some rational velocities in (ii,si) and the other 
is the structure of P in (i 2 ,s 2 ) \ when the set of increase points of h 2 is 

bounded. 

In order to have these both characteristics, we set for a G (0, 1) and e G (0, 1) 

L 1 = ^(a 1 +a- 1 ) and L 2 = (1 - a) ((1 - e)Id + | {a 2 + a" 2 )) , (21) 

where Id is the identity on M. It is immediate to verify that these obey the 
linearity continuity, positivity, s-homogeneity and normalization conditions. It is 
also easy to see that D(h x ) = E(hx) = {-1, 1} and D(h 2 ) = E{h 2 ) = {-2, 0, 2} and 
consequently E{h 2 ) DR(E(hi)) = 0. Therefore, Proposition |3j] does not apply and 
some rational velocities are expected to be absent from P. The results are contained 
in the following statement. In the remaining of this section, p and q denote two 
co-prime integers. 

Proposition 5.1 For the couplings given by (21) we have P C [—2,2] and 

(i) the velocity v G [— 1, 1] belongs to P iff either v is irrational or v = 2 with p + q 
odd, 

(ii) the velocity v G [1,2] (resp. [— 2, — 1]) belongs to P iff it exists k G N such that 
v = 1 + FFI ( res P- v = -l- 

Notice that the set 

f P 

[— \ P = < v £ [— 1,1] : v = —, p + q even 

I 1 

is dense in [— 1, 1] but G is nowhere dense. 

We have plotted on Figure 2 the selected velocity v(c) versus c G [|, 1], the func- 
tion being symmetric, i.e. v(c) — —v(l — c + 0). This function is a Devil's staircase 
in [v~ (— l),w _1 (l)] and a classical staircase in the complementary intervals. 

Proof of Proposition 5.1: (i) By a direct calculation, we obtain 

00 

<j>(x, V) = (1 - a) Yl (( J " e )'« + + H n-v(k+l) (x) Vx G M, 

fc=0 nez 

where the coefficients 1* are the entries of the powers L\. These coefficients are 
known explicitly . In particular they obey the following property 

l\ ± iff |n| < k and G Z. (22) 



On one hand all the irrational velocities [—1,1] are possible as given by Lemma |3.3| 
(i). On the other hand we have 



(a;, = 0(0 - 0, -) = 0(0, - - 0) for x G 

9 q q 



i.o 

9 



and thus the proof of Theorem ^2| tells that | € [-1, 1] belongs to P iff 0(0, | -0) < 
0(0, |). Using the definition of the Heaviside function, we have 

oc 

0(0, H) - 0(0, £ - 0) = (1 - a) £ ((1 - e)^- 1 + |(C=2 + 

11 k=1 K L ' 
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Figure 2: The selected velocity v(c) versus c for the one-dimensional lattice dynam- 
ical system example (a — e = 0.6). 



Taking into account relation (|22|), it turns out that 

0(0,^) = 0(0,^-0) iff Vfc > 1, 3j £ N : k(p + q) = 2j, 
from which we conclude that | £ P iff p + q is odd. 



(ii) is easily deduced from the proof of statement (ii) Lemma 3.3 



□ 



5.2 Continuous example 

As a continuous couplings example, we consider for a £ (0, 1), L\ = aL and L 2 = 
(1 — a)L where L is now the continuous diffusion operator expressed in the integral 
formulation, i.e. the corresponding function h is defined by 

h{x) = ^ e-^dy, 

so that the couplings obey the linearity, continuity, positivity, s-homogeneity and 
normalization conditions. In this case we have E{h\) = E(h2) = K and D(hi) = 
DQ12) = thus P = M and G is empty. We conclude that for any a,c £ (0, 1) it 
exists a unique front shape. The corresponding velocity is a continuous function of 
(a,c) with range M. 

The front velocity v(c) in this continuous model, restricted to the values in [0,1], 
is shown on Figure 3. Notice that v(.) is an odd function. 



5.3 Multi-dimensional Coupled Map Lattices 

We now turn to the application of the previous results to the dynamics of pla- 
nar fronts in discrete bistable extended systems defined on a d-dimcnsional lattice, 
by considering a simple situation, that of homogeneous CML. However a similar 
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Figure 3: The selected velocity v(c) versus c for the continuous coupling (a = 0.4). 



analysis can be applied, without supplementary difficulties, to more general multi- 
dimensional extended mappings. After defining these systems, we shall see how the 
homogeneity implies the reduction of the dynamics, for the orbits having equal com- 
ponents on (d— l)-dimensional planes, to a one-dimensional extended mapping. The 
existence of planar fronts in these multi-dimensional models then directly follows 
from this reduction and from Theorem 3.2. Finally, the monotonicity of the dynam- 



ics allows to conclude about the velocity of interfaces, even when their dynamics 
does not reduce to a one-dimensional extended mapping. 



5.3.1 Definitions 

Multi-dimensional CML are introduced following the ideas of |J. For d G N, d > 0, 
let 

M' = {u = {u s } seZ d, it s e R : \\uWoo < +00} , 

be the phase space under consideration. A homogeneous e?-dimensional CML is a 
discrete time dynamical system for which each time iteration of an orbit {it*} of 
elements of M! is given by the following relation 

u t+1 =LFou t . (23) 

Here the mapping F is a direct product of one-dimensional identical mappings, i.e. 

{Fou) s = f(u s ) VseZ d , 

where / maps R into itself. The coupling L is required to be a linear, continuous, 
positive and s- homogeneous operator of M! into itself. Its explicit expression is 
then 

(Lu) s = ^ l n u s - n Vs G Z d , 
where the real coefficients l n ^ obey the normalization condition 
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A point s in the multi-dimensional lattice 7L d is denoted by {siji^^d- The dif- 
ference of two elements in 7L d is to be understood as the component-wise difference. 
Moreover we are going to use the Euclidean inner product in R d , i.e. 

d 

r.s = TiSi Vr, s € R d . 



5.3.2 One-dimensional orbits 

Multi-dimensional homogeneous CML satisfy any translational invariance symmetry 
along each lattice axis. To each of these symmetries corresponds a reduction of 
the dimension of the lattice supporting the dynamics. Hence by combining d — 1 
linearly independent translations, the resulting orbit's dynamics is governed by a 
one-dimensional extended mapping, as we describe now. 

Let k £ R d be such that k.k = 1, a fixed direction. By induction and by using 
the coupling homogeneity, one shows that for any orbit {w* }tgN generated by the 



dynamics (23), if the property 

Vr, s E Z d : k.r = k.s u\ = u\ , (24) 

holds for t = 0, then it is satisfied for any t 0. In this case0 the time evolution 
reduces to that of a one-dimensional discrete extended mapping. Actually by intro- 
ducing the reduced coordinate u = k.s and the corresponding component z l u = u\, 
we obtain the component-wise expression for the reduced dynamics 



z t+1 



E VcoeZ(fc), (25) 



where 

Z{k) = {uj e K : u = k.s, s 6 1 d ) and A w = ^ /„ Vw 6 Z(k). 



Similarly as in Remark 3J3, if Vz , j £ Q then it exists m £ R + such that Z(fc) = 



mZ and if Eli, j such that p £ K \ Q then is dense in R. Furthermore for a 

given multi-dimensional coupling L, the one-dimensional (discrete) reduced coupling 
A defined by 

uez(k) 

and thus the dynamical properties of the corresponding one-dimensional orbits of 



(23), depend on the direction k. From now on we choose as a local map, the piece- 
wise affine model (j|) . The reduced system (^5|) then turns out to be a particular case 
of the system (Q) with the coupling functions given by hi — ah k and hi = (1 — a)h k 
where 

h k (uj) = ^2 KH (lu - v) = l n H Q (uj - k.n). (26) 

u£Z(k) nez d 

5.3.3 The planar fronts 

Definition 5.2 A planar front of velocity v in the direction k is an orbit of 1^2- 
defined by 

u* = <j>(k.s - vt - x ) VteZ,seZ d , 
where xq £ R and the shape 4> £ M. and obeys the condition ^). 

10 This restriction is not necessary if k is an irrational direction, i.e. if k.s = implies s = 0. 
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As for fronts, the existence of planar fronts only depend on the existence of their 
shape, and is therefore independent of the translation parameter Xq. 

Since the planar fronts are invariant under translations orthogonal to their di- 
rection, we can apply the previous reduction to obtain the equivalence between 
the planar fronts of ( |2^ ) and the fronts of the corresponding reduced system ( p5| ) . 
The planar front existence in piece-wise affine bistable d-dimensional CML then 
follows from Theorem 3.2 and the propositions of Section 0. Let us denote by 



i k = inf {uj e Z(k) : A w > 0} and by s k = sup {lu e Z(k) : A w > 0}. 

Corollary 5.3 Given a, L and k, there is a countable nowhere dense set G such 
that, for any c £ (0, 1] \ G, there exist planar fronts in the system (21) with a unique 
shape. The corresponding velocity is a continuous function of (a,c,k,L) with range 
[ik, Sk]- 

In particular for a given coupling, the selected velocity as well as the extreme 
velocities and the plateaus, depend continuously on the planar fronts' orientation. 

As a sketch of proof, one firstly notices that from ( p5| ) it follows the relations 
E(fa) = D(hx) = E(h 2 ) = D(h 2 ) C Z{k) n [i k ,s k ] and hence P = E(<p(0,.)) = 
[ik, s k ] and G — G 2 - This allows to state the existence following the existence's proof 
in the general case. The continuity versus k and L is a consequence of Proposition 



3.11 and of the fact that small variations of k and of L induce small variations of 



the reduced coupling A for the distance d(., .) as can be easily seen from fl26|). 
5.3.4 Velocity of interfaces 

Similarly as in the one-dimensional case, the subsequent question to the existence 
of fronts is the velocity of interfaces in multi-dimensional CML. Naturally their 
definition must take into account the orientation. 

Definition 5.4 A (fc, c) -interface is a configuration u S M! such that 
3J ki ,J ks eR -[ Us ^ c if sk>Jks ■ 

The set of (fc, c) -interfaces is denoted by J k and 3 k \{u) and J ks (u) are the functions 
from J k to K defined by 

Jki{u) = inf{s.k : u s ^ c} and J ks (u) — supjs.fc : u s < c}. 

To each (fc, c)-interface u we can associate two (fc, c)-interfaces u and u satisfying 
relation (|2~i| ) and defined for any s 6 Z d by 

u s = inf{u„ : n.k — s.k} and u s — sup{u„ : n.k = s.k}. 

We have u s ^ u s ^ u s Ms 6 Z d and then J k i(u) ^ Jki(u) ^ Jfei(u) and J ks (u) ^ 
Jks{u) ^ Jks{u)- For the configurations u and u the dynamics can be reduced to a 
one-dimensional extended mapping as was done for planar fronts. Consequently we 



can apply Proposition 4.2, Theorem 4.3 and the order preservation of the dynamics 
to conclude that a (fc, c)-interface initial condition vP generates a (fc, c)-interface 
orbit {u*}t<=N and the following statement holds. 

Corollary 5.5 For any u° € J k c the following limits exist and do not depend on 
the particular choice of u° 

Jkiiu*) Jksiu*) 
lim — lim = Vk, 

t^+oo t t^+oo t 

where Vk is the velocity of the planar front in the direction k. 
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6 Other fronts and general bistable extended map- 
pings 

Starting from the knowledge of the fronts dynamics in piece-wise afhne bistable 
extended mappings with homogeneous linear couplings, one can extend the results 
to other kinds of fronts and to more general systems in which one or many of the 
previous conditions are relaxed. Here we firstly consider changes in the local map 
and then we allow for nonlinear couplings. 

Let us recall that a fronts velocity v corresponds to a plateau iff it is a dis- 
continuity point of the function 0(0, .). In this case according to jis|), the interval 
corresponding to this plateau is given by /.„ = (0(0 — 0,v), 0(0, v)], the left bound 
being included if v £ Pi- In the following I® denotes the interior of this interval. 

In the interior of a plateau and without changing the couplings, the existence of 
fronts is claimed in systems defined by (|l|) but with now a C°° bistable local map. 
Actually inside the plateaus, the distance from the front shape to the discontinuity 
c is positive. In this case one can modify the local map in a neighborhood V of 
c where the fronts never go, i.e. V C Zj, without modifying these orbits. The 
existence of fronts in some bistable extended mappings with a C°° local map then 
directly follows from this observation. Similarly it is possible to extend the result 
on the velocity of interfaces to these mappings. 

Condition (0) in the front definition imposes the shape to be a right continuous 



function as it is claimed in relation (11). However one can ask about the existence 
of fronts of a different kind, namely fronts having a left continuous shape satisfying 
instead of (||), the condition 

l(x) < c if x < 
l(x) ^ c if x > 

Obviously, if such fronts exist, we have 

0(0) = <f>(x - 0), 



where is given by (10). Consequently according to Theorem 3.2, the dynamical 
system (^) admits fronts with a left continuous shape iff it has fronts with a right 
continuous shape which does not touch c (i.e. 0(0) ^ c). The latter occurs when 
the velocity belongs to the interior of a plateau. Moreover following the previous 
argument, fronts with a left continuous shape also exist in a larger class of bistable 
extended mappings for which smooth local maps are allowed. 

Still in J° for v £ D(0(O, .)), the local map can be modified (now only) at c in 
order to obtain a third kind of fronts still with the same velocity v but which will be 
unstable (in the sense of Lyapunov) .0 One has to consider the following definition 
of the Heaviside function instead of the one used up to now. 

( if x<0 
H r Q {x) = \ r if x = 
{ 1 if x>0 

for Q ^ r ^ 1. The corresponding front shape follows from the introduction of 
this function in condition (^|) and from solving the subsequent shape equation. 
The solution now depends continuously on r and we have, writing explicitly the 
dependence on r, 

4>'{0 - 0, v, 1) = </>'(0, v, 0) < </>'(0, v, r) < (f>'(0, v, 1) = 0(0, v). 



11 Notice that Theorem 
velocity. 



4.3 imposes all the fronts defined here to propagate with the same 
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Hence given c inside a plateau, the front shape <fi' exists if r is such that 

4>'(0, v , r) = c and /(c) = ac + r(l — a). 

It is unstable since any small enough positive (resp. negative) perturbation of at 
belongs to the basin of attraction of the front with shape cf> (resp. <j>). Actually 
one can show the linear stability of the latter fronts inside the plateaus. For r = 1, 
4>' coincide with (f> whereas for r = it is nothing but 4>. Moreover the extension 
from discontinuous to C°° local map can still be done inside a plateau by requiring 
/(c) — ac + r(l — a). The front shape built with the present definition of the 
Heaviside function will be stable if |/'(c)| < 1 and unstable if |/'(c)| > 1. 

According to these results, we conjecture that in the dynamical system (|l|), if 
the local map is continuous and bistable, there are fronts with right continuous 
shape, fronts with left continuous shape and unstable fronts, the various shapes 
only differing at their discontinuity points. 

Keeping on the assumption c 6 /J, one can modify the couplings in the following 
way. We perturb C and £ with two homogeneous operators in M, namely C and 
£, of class C 1 . Instead of ([]]) we now consider the following mapping on M. 

u t+1 =T{u\e), (27) 

where 

T{u, e) = ((1 - e)C + e£ju+ ((1 - e)£ + e£') F o u. 

A front of velocity v in this dynamical system is a travelling wave for which the 
shape obeys the equation 

Q{u, e) = F(u, e) — cr v u = 

with the condition H c ou = Hq (where the Heaviside function is the right continuous 
one). Its existence can be proved using a continuation of the shape ^> at e = 
(similarly as for the continuation of solutions from the uncoupled limit in p"T|). 

Proposition 6.1 Given C, C! , C and £ , v £ D(<p(0, .)) and c 6 7°, it exists €q such 
that for all e £ [0, e ), the system has fronts of velocity v; their shape being a 
continuous function of e. 

Proof: For e = let <f> be given by Theorem |3.2| . By assumption on v and c we have 

inf \4>(x) — c| = S, 

X 

where S is given by 

S = min {c - 0(0 - 0, v), 0(0, v)-c}> 0. 
Let U be the neighbourhood of <\> in M. defined as the set of functions u such that 

3x £ K : \\</)-a xo u\\ < S, 

hence for u £ U, we have 

J u{x) < c if x < x° 
\ u(x) > c if x ^ x° 

In U, Q is C 1 since C, £ and F are. Moreover D u Q(cf),0) is invertible and its 
inverse is bounded as follows from solving the front shape equation in (|l|). By the 
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implicit function theorem it exists eo and a unique continuous function u(e) such 
that u(0) = 4> and for any e G [0, eo), we have it(e) £ U and Q(u(e), e) = 0. □ 

Finally we notice that no supplementary difficulties are found when one combines 
the two previous extensions, that is to say when one both perturbs the local map 
and the couplings. 
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A Proof of Lemma 13. 3 



In all the proof, we will employ the sequences {fcnjneN, k n € N and {e n }neth £n S K, 
e n > so that 

lim k n = +oo and lim e„ = 0. 

n — >+oo n — >+oo 

(i) By assumption let a,b E E(hi) be such that a < b and let d £ E(li2)- Since is 
irrational the sequences fc„ and e n can be chosen so that 

-e„<- j(fc„ + l)6»-^— ^| SCO VneN, 

where {.} stands for the fractional part of a real number 0. Moreover since < 
9 < 1 the sequence 

j n = L(fc n + i)0-^j, 

satisfies the inequalities ^ j„ *S &n for ft sufficiently large, say n ^ Uq. Here 
denotes for the floor function [Q. Using the equality x = \x\ + {x}, we obtain 

-(b - a)e n < d+ (k n - j n )a + j n b - (k n + 1) (a + 6(b - a)) < Q Vn ^ n , 

that is to say a + 9(b — a) G P. 

(ii) From relation (nj) let us assume that 



^ fc + 1 

fc=0 



Then there exist the sequences {a^njngN.i^i^fcn , a i,n G E(hi), and {e^l^gN, <^n G 
E(h 2 ) such that 



-Cn < , 3- : ( + 0>i,n ) - V 

fcn + i V i=i / 

Using the definition of and Sj, we have 



- v < VneN. 



Consequently if E(h\) and EQ12) are bounded then i x ^ u ^ &\- 111 this case one 
has 

Em .)) n Sl r = u n [ h , SlT , 

and fiij is proved. 

(Hi) For > 1, 6 irrational and v = a + 6(b — a) where a, b <E E(h\) are such that 
a < b, let {<i n }„gN be an increasing sequence of elements of E(h<z) such that 

lim d n = +00. 

n — >+oo 

Let us define 

Pn = -r and k n = \-^\ , 

— a y 

where \x] stands for the ceiling function |?]|. By compactness it exists a £ [0,9] 
and a strictly increasing sequence of integers {rii}igN such that 

a = lim k 7li 9 — (3 ni . 

i — >+oo 



2G 



Let also {rrij}j0$, frij G N be such that 

j<{mjfl + a}<-^. 
The increase of /3„ results in the following. For all j it exists such that 

-1 + -f— x f - > mj and \k Ui . - /3„ 4 . - a| < -± . 

Finally let kj = k ni + rrij . From the previous inequalities, we obtain the following 

sc lkj6-p ni .\ < kj 

and 

< {kjO-pn,.} < €j. 

Now let pj — [kjd — (3 ni J to obtain 

-(& - a)€j < dm. + (kj - pj)a+pjb - (kj + 1) (a + 6(b - a)) ^ 0, 
and then to conclude that v G P. 

B Proof of Lemma |5.1U| 

Lemma 3. 10| is a consequence of the combination of the following statements. 



Lemma B.l (i) If D(Ji2) ^ and E(hi) = D(h\) then L>(</>(0, .)) is dense in 

(«l,Sl). 

(ii) If sup D(li2) = +oo and a G D(h\) then D(cf>(0,.)) is dense in (a,+oo). 
(Hi) If inf D(h2) = — oo and a G D(h\) then D((f>(0, .)) is dense in (— oo, a). 

Lemma B.2 If D(cf)(0, .)) is dense in (i>i,t>2) then G is nowhere dense in v~ 1 ((vi,V2))- 
Proof of Lemma |B.i[- Let a,b £ D(hi) and d G D(h2)- We have 



vo<;<* d+(fc - + i) 1 a+i6 6^(o,.))- 

Hence for any G [0, 1] by choosing j/j = [fc#J we obtain 

d+ (fc -jk)a + jkb 



a + 8(b — a) = lim 



fc + 1 



which proves (i). 

Now let {d n } ne N, d n G D(h,2) and lim (i„ = +oo. For 6* > by choosing k n = I 4f J 

n — >oo 

we have 

%±^£ e 15(^(0, .)), 

and thus 

a + u = urn — — , 

k— >oo ft n -|- 1 

(nj is proved and one can proceed similarly for (Hi). □ 



To prove Lemma B.2 ( we need the following auxiliary lemma. 
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Lemma B.3 If v{c\) — v{c2) then (01,02) nG = 0. 



Proof: Let c\ < c 2 be such that v(c\) = vfa). For any c £ (ci,C2) we obtain, by 
using the definition of the function v, v(c) = v(c\) = v(c 2 ) and 

0(0, v(c)) > c 2 > c> ci ^ 0(0, u(c) - 0) = 0(0 - 0, w(c)), 

and hence c ^ G. □ 

Now let ci < C2 be such that ci,C2 6 G fl u _1 ((wi, «2))- Then either v(ci) — 
v(c 2 ) and by Lemma |B.3| (ci,c 2 ) D G = 0, or v(ci) < v(c 2 ). In the latter case 
by density of Z?(0(O, .)) it exists v* e (v(a),v(c 2 )) fl Z?(0(O, .)). Let us choose 
c*,^ e (0(0, u* - 0), 0(0,1)*)) s uch that c\ < c 2 . We have (c*,^) C (ci,c 2 ) and 



w(c^) = ^(cj). So again Lemma B.3 implies that (c*, c*.) (1 G = 0. In both cases G 



is not dense in (c\,c%). □ 
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